Cosmology of a covariant Galileon field 



O 
(N 

< 



o 
u 

Oh' 
I 

o 



(N 
> 

O 

o 

(N 

l> 

O 

o 



% 



Antonio De Felice and Shinji Tsujikawa 

Department of Physics, Faculty of Science, Tokyo University of Science, 

1-3, Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan 

(Dated: August 13, 2010) 

We study the cosmology of a covariant scalar field respecting a Galilean symmetry in flat space- 
time. We show the existence of a tracker solution that finally approaches a de Sitter fixed point 
responsible for cosmic acceleration today. The viable region of model parameters is clarified by 
deriving conditions under which ghosts and Laplacian instabilities of scalar and tensor perturbations 
are absent. The field equation of state exhibits a peculiar phantom-like behavior along the tracker, 
which allows a possibility to observationally distinguish the Galileon gravity from the ACDM model. 



The problem of dark energy responsible for cosmic ac- 
celeration today has motivated the idea that the grav- 
itational law may be modified from General Relativity 
(GR) on large scales (sec Refs. [l[ for reviews). On the 
other hand, one needs to recover Newton gravity at short 
distances for the compatibility with solar-system exper- 
iments. Besides the chameleon mechanism [2| based on 
the density-dependent matter coupling with a scalar field 
(used also in f{R) theories 0,3), there is another way 
to recover GR at short distances: the Vainshtcin mech- 
anism |5| based on non-linear field self-interactions such 
as n0(V0)^, where (V0)^ = d'^4'dfj,4>. This non-linear ef- 
fect has been employed for the brane-bending mode of the 
self-accelerating branch in the Dvali-Gabadadze-Porrati 
(DGP) braneworld [y], but the DGP model is unfortu- 
nately plagued by a ghost problem Q. 

In order to avoid the appearance of ghosts, it is im- 
portant to keep the field equations up to second-order in 
time-derivatives. A scalar field called "Galileon" Q, 
whose action is invariant under the Galilean symmetry 
df^cj) — ?> dfj.(j) -\- bf^ in flat space-time, allows five field La- 
grangians that give rise to derivatives up to second-order 
(see Refs. [a|-[ll| for related works). If the analysis in 
[3| is extended to the curved space-time, one needs to 
introduce couplings between the field (f) and the curva- 
ture tensors for constructing the Lagrangians free from 
higher-order derivatives in the equations of motion Q . 

The five covariant Lagrangians that respect the 
Galilean symmetry in flat space-time are given by 



£1 = M^ 



C2 = (V0)^ 

\2 



£4 = (V0)2 [2{acj)f - 2cj),^A-'"' - i?(V(/.)V2] /A/^ 
£5 - (V0)2[(n0)3 - 3(00) 0;^.0'^'' 

+ 20;/(/);,''0;/ ~ 6 ^ ; ^ ^ ' ^"^ ^ ^ ^ G,p ] / Af ^ , (l) 



where a semicolon represents a covariant derivative, M 
is a constant having a dimension of mass, and Gi^p is the 
Einstein tensor. In this Letter we study the cosmology 
based on the action 



f^^h^^l 



d X Cm . 



(2) 



where Mpi is the reduced Planck mass and Ci are con- 
stants. For the matter Lagrangian Cm we take into ac- 



count perfect fluids of radiation (density pr) and non- 
relativistic matter (density pm)- Although the cosmolog- 
ical dynamics up to £4 were discussed in Ref. (11| . we 
will show that inclusion of £5 is crucially important to 
determine the full Galileon dynamics. Moreover the vi- 
able parameter space will be clarified for such full theory. 
In the flat Friedmann-Lemaitre-Robertson- Walker 
(FLRW) Universe with a scale factor a{t), the variation 
of the action (U) leads to the following equations of mo- 
tion 






'pH 



^pi^ 



-Pde 



Pr/3, 



(3) 
(4) 



where H = a/ a is the Hubble parameter (a dot represents 
a derivative with respect to cosmic time t), and 



PDE = -ciA'/3^/2 - C20V2 + 3c3i/0VA'/^ 
-45c4i?20V(2MS) + 21c5H^(j>yM^, 
'de = ciM30/2 - C20V2 - C30'0/A/3 

+3cA^^[8H^ + (3^2 ^ 2H),j)]/{2M^) 
-3c5H^^[5H4> + 2{H^ + H)<j>]/M^ . 
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(5) 



(6) 



The matter fluids obey the continuity equations p„j 4- 
3Hpm = and pr + ^Hpr = 0. From Eqs. (O and (g]) the 
dark component also satisfies /Ode + iH{pY)E + -Pde) = 0. 
We define the dark energy equation of state wde and 
the effective equation of state WcS, as wde = Pde/ Pde 
and Wcff = — 1 — 2H/{3H^). The latter is known by the 
background expansion history of the Universe. 

Since we are interested in the case where the late-time 
cosmic acceleration is driven by field kinetic terms with- 
out a potential, we set ci = in the following discussion. 
When ci = and C2 ^ the only solution to Eqs. ([3]) 
and (U]) in the Minkowski background {H = 0) without 
matter corresponds to = 0. We introduce the following 
quantities useful to describe the cosmological dynamics 



n = <j>dsHds/m) , 



r2 



(0/0ds)Vn 



(7) 



where 0ds and H^s ~ 10~^° Mp\ are the field velocity 
and the Hubble parameter at the de Sitter (dS) solution, 
respectively. The mass M is related to H^s via AI^ = 
Afpii^lg. At the dS point one has ri = 1 and r2 = 1. 



Equation (j3]) can be written as Vim + ^r + ^de = 1, 
where n^ = p^/iSM^^H^), n, = p^ / {'iMl^H^) , and 
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^DE = --C2x\Qrlr2+c-ix\srlr2 — —Cix\^rir2+lc^x%^r2, 

(8) 

where Xds = '?^ds/(-ffds-^^pi)- Since J^de = 1 at the dS 
— 2 

dS' 



point, Eq. ([5]) gives a relation between the terms C2X 



caSjg, a = C4x|^g, j3 = csx^g. Combining this with an- 
other relation coming from Eq. Q, we obtain 

C2xl^ = 6 + 9a - 12/? , caxi^s = 2 + 9a - 9/3 . (9) 

It is useful to use a and /3 because the coefficients of 
physical quantities (such as flDE) can be expressed in 
terms of those quantities thanks to Eq. ^ . The relations 
^ are not subject to change under a rescaling Ci — )- Ci/7* 
and Xds -^ IXds, where 7 is a real number. Therefore 
rescaled choices of Ci will lead to the same dynamics (as 
they have the same a and P) for both the background and 
the linear perturbation, which implies that redefining the 
coefficients Ci in terms of a and (i is convenient. 

The autonomous equations for the variables ri , r2 , fir 
follow from Eqs. (O, (JH), and fluid equations. One can 
show that there is an equilibrium point characterized by 



ri = 1 , i.e. 4>H = constant , 
at which the variables r2 and fir satisfy 



2r2 (3 - 3r2 + »r) 
l + r2 



n:.^ 



nr {Qr-l- 7r2) 

l + r-2 ' 



(10) 



(11) 



where a prime represents a derivative with respect to 
N = In a. This result is interesting because it shows 
the universality of the equations of motion without any 
dependence on a and /3. Along the solution ((TU)) . the field 
velocity evolves as ex t during radiation and matter 
eras {H oc 1/i). There is also a simple relation 51de = 7*2 
along the solution n = 1. 

We have three fixed points: (a) {ri,r2,flr) = (1,0,1) 
[radiation], (b) {ri,r2,i^r) = (1,0,0) [matter], (c) 
(ri,r2,rir) ~ (171,0) [dS]. The stability of these points 
can be analyzed by considering linear perturbations 
6ri , dr2 , SQr about them. The perturbation Sri satisfies 



6r[ = - 



9 + »,. + 3r2 
2(1 + r2) 



6ri , 



(12) 



which shows that, in the regime < r2 < 1 and O^ > 0, 
the solution is stable in the direction of ri . Since the dS 
point is stable in the other two directions, the solutions 
finally approach it. The points (a) and (b) are saddle 
because they are unstable in the direction of r2. 

Along the solution dTU]) we have pde = 3M'^/H^, 
Pde = -3Af6(2 + w,s)/H^, and 



w-DE = -2-ti;off = - 






3(r2 + l)' 



WcS 



D,r — 6r2 
3(r2 + l)' 



(13) 



From the radiation era to the dS epoch the effective equa- 
tion of state evolves as WeS = 1/3— >0— >■— 1, whereas 



the dark energy equation of state exhibits a peculiar evo- 
lution: wde = -7/3 -> -2 -> -1. 

The evolution of dark energy is different depending on 
the initial conditions of (ri, r2, fir). If they arc chosen to 
be close to the fixed point (a) at the onset of the radiation 
era, then the solutions follow the sequence (a) — > (b) — > 
(c). If ri ^ 1 initially, the dominant contribution to 
r^DE comes from the term £5, i.e. fluE — 7/3r2. In this 
case the solutions approach ri = 1 at late times with the 
increase of ri . For the initial conditions with ri 3> 1 the 
term £2 gives the dominant contribution to J^de, but 
this case is not viable because the field kinetic energy 
decreases rapidly as in quintessence without a potential. 
Numerical simulations show that if ri < 2 initially the 
solutions approach ri = 1, but in the opposite case the 
Universe finally reaches the matter-dominated epoch. 

Let us find the allowed parameter space in terms of 
(a, /3) by deriving the conditions for the avoidance of 
ghosts and instabilities of scalar and tens or p erturba- 
tions. Using the Faddeev-Jackiw method [l2l, the ac- 
tion (12) can be expanded at second-order in perturba- 
tions. Following the similar procedure to that given in 
Refs. [1^ , the no-ghost condition for the scalar sector of 
the action is given by 



Qs = -s/(l-fM3)'>0, 



(14) 



where s — 6(1 + /Ui)(/Ui + P2 + Mi/^2 - 2/i3 — pi), and 

Ml = 3arir2/2 - 3/3r2 , (15) 

P2 = (3a ~Af3 + 2)rlr2/2 - 2(9a - 9/3 + 2)r2r2 

-I- 45arir2/2 - 28;3r2 , (16) 

P3 = -(9a - 9/3 -h 2)r2r2/2 + 15arir2/2 - 21;3r2/2. (17) 

The condition for the avoidance of Laplacian instabilities 
associated with the scalar field propagation speed is 

4 = {(1 + pi)^[2p'3 - (1 + P3){5 + 3wcff) + 3f^„, + 4f}^] 
-Ap[{l + Mi)(l + P2) + 2(1 + P3)^il + Pi)}/s > 0, (18) 

where 

Pi = -arir2/2+3/3r2(3+3u;cff-3r'i/ri-r^/r2)/2. (19) 

Similar calculations for the tensor perturbation lead to 



Qt = 3arir2/4 - 3^r2/2 + 1/2 > , 



(20) 



2 ^ 2ri (2 - arir2) - 3/3 {r2r[ + nr'^) , 

^ 2ri(2-f 3arir2-6/3r2) '^ ' 

We consider the following three different regimes. 

• (i) ri < 1, r2 < 1 

This characterizes the early cosmological epoch in 
which the term C5 dominates the dynamics of the 
field. For the scalar modes we have Qs — 6Q/3r2 
and c| ~ (l-|-r2r)/40. The sign change of r2 implies 
the appearance of ghosts. For the initial conditions 



J a= 12/3"--9i3-2, cr- = Oat(cr-)' = 



a = 2/3+-, (Gs),,,, = 







;8 



FIG. 1: The viable parameter space in the {a, /3) plane for the 
branch r2 > 0. We also show several conditions that deter- 
mine the border between the allowed and excluded regions. 



with r2 > 0, it is required that /3 > 0. The Lapla- 
cian instabilities of the scalar modes can be avoided 
because c| ~ 1/20 and c| ~ 1/40 during radiation 
and matter eras, respectively. Since Qt — 1/2 and 
c^ ~ l + 3(3r2{5 — 3Qr)/& — 1, the tensor modes do 
not provide additional constraints. We also have 

K;DE^-(l + a-)/8, Wcff=if^r/3, (22) 

which is valid for fir 3> {ri, r2\. 

• (ii) ri = 1, r2 < 1 

This corresponds to the equilibrium point (|10p dur- 
ing radiation or matter domination. The conditions 
(0 and dH]) reduce to 



Os ^ 3(2 - 3a + 6/?)r2 > , 
2 _ 8 + 10a - 9/? -I- VLr{2 -I- 3a - 3/3) 
^^~ 3(2 - 3a -h 6/3) 



(23) 
>G. (24) 



For the branch r2 > the first condition reduces 
to 2 — 3a -|- 6/3 > 0. For the tensor modes, we have 
4 ~ 1 - r2(4a -h 3/3 -H 3/3f7^)/2 ~ 1 and Qt > 0. 



1, r, 



1 



• (iii) ri 

This corresponds to the dS point, at which the con- 
ditions (HH), (gni), (0, and dni) are given by 



_ 4-9(a-2/3)^ 
^^- 3(a-2/3)2 >0' 
Qt = (2 + 3a - 6/3)/4 > , 
2 _ (a - 2/3)(4 -}- 15a2 - 48a;3 + 36^^) 



C'-p 



2-a 



2 -f 3a - 6/3 



2[4-9(a- 2/3)2] 
>0. 



(25) 

(26) 

> 0,(27) 

(28) 



If /3 > 0, c|n can have a minimum during the transition 
from the regime r2 ^ 1 to r2 ~ 1. This value tends to 
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FIG. 2: Evolution of WcB and iude for the cases: (A) a = 
— 1.4, (5 = —0.8, Xds = 1 with initial conditions r\ = 1, 
r2 = 10"'*°, fir = 0.99999 at the redshift z = 3.11 x 10^ 
(B) a = 0.1, /? = 0.049, Xds = 1 with initial conditions r\ = 
5x10"", r2 = 8xl0"^^!:^r =0.999995 at z = 6.44xl0^and 
(B') the same a, /3, idS as in the case (B) but with different 
initial conditions ri = 5 x 10"^, r2 = 8 x 10"^'', Q.r = 0.9995 
at 2 = 6.72 X 10^ 



decrease as /3 approaches 1. Imposing that c|n > at 
the minimum, we obtain a < 12v^ — 9/3 — 2. In Fig. [T] 
we illustrate the region in which this condition as well as 
the conditions ((23)) -p8 l) are satisfied for r2 > 0. Numer- 
ical simulations confirm that for the parameters inside 
the shaded region in Fig. [T] the no-ghost and stability 
conditions are not violated even in the intermediate cos- 
mological epoch. 

In Fig. [2] we plot the variation of wde and Wcs versus 
the redshift z for several different model parameters and 
initial conditions. In the case (A) the initial condition is 
chosen to be ri = 1, so that wde and Wcs evolve accord- 
ing to Eq. ((T5)) with the variation of 7*2 and fir ■ While the 
evolution of WeS is similar to that in the ACDM model, 
the dark energy equation of state evolves from the regime 
wde < — 1 to the dS attractor with wde = ^1- The cases 
(B) and (B') in Fig.[2]correspond to the initial conditions 
in the regime (i). As estimated by Eq. (P^ . wde starts 
to evolve from —1/4 and reaches the value —1/8 during 
the matter era. The evolution of wde is different de- 
pending on the epoch at which ri grows to the order of 
1. In the case (B) the solutions reach the regime ri ^ 1 
only recently, whereas in the case (B') the approach of 
this regime occurs much earlier. The equilibrium point 
PU)) can be regarded as a tracker that attracts solutions 
with different initial conditions to a common trajectory. 
Before approaching the tracker, the solutions cross the 
boundary wde = — 1 without any unstable behavior of 
perturbations. Note that there are no significant differ- 
ences for the variation of Wcff between the cases (A) and 
(B) [also (B')]. 

In Fig. [3] we show the evolution of Cg and c^ for the 
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FIG. 3: Evolution of cg and Cy for the cases (A) and (B) as 
in Fig. [2] 



same model parameters and initial conditions as in Fig. [5] 
In the case (A) the scalar propagation speed remains sub- 
luminal, as estimated by Eqs. (P^ and ((?7)) . For a ~ 
-1.4 and /? = -0.8, Eq. ^Q shows that at the dS point 
the tensor mode becomes supcr-luminal. However, both 
the scalar and tensor modes can be sub-luminal at the dS 
point, as in the case (B) of Fig. [3] (a = 0.1, (3 = 0.049). 
For the initial conditions starting from the regime (i) 
we require /3 > to avoid ghosts. Under the conditions 
p3)) . ((25|) . and ^ with /3 > 0, one can show that c| in 
Eq. (P^ gets larger than 1. If the solutions approach the 
tracker in the regime (ii) long before the dS epoch, there 
is a period in which c| exceeds 1. This super-luminal 
propagation can be avoided if ri grows to the order of 
unity only recently. The case (B) in Fig. [3] corresponds 
to such an example for which c|. has a peak smaller than 
1 after the matter era. In this case the tensor mode 



is slightly super-luminal in the regime (i). In general 
there is a period in which the propagation speed of either 
scalar or tensor modes exceeds 1. However, this does not 
necessarily imply the inconsistency of theory because of 
the possibility for the absence of closed causal curves [Il[ . 

In summary we have studied the cosmology for the full 
Galileon action ([2]) and derived all conditions for the con- 
sistency of such theory. We have shown that, under the 
conditions ([9]) , there exist stable dS solutions responsible 
for dark energy. In spite of the complexities of Galileon 
Lagrangians, the conditions for the avoidance of ghosts 
and Laplacian instabilities constrain the allowed param- 
eter space in terms of the variables a and /3 in a simple 
way. While the evolution of iude, <?s ^^"^ ^t i^ different 
depending on the model parameters and the initial condi- 
tions of ri , we have derived convenient analytic formulas 
to evaluate those quantities in three distinct regimes. 

There are several interesting applications and general- 
izations of Galileon gravity. First, the study of cosmo- 
logical perturbations may provide some signatures for the 
modification of gravity from GR. The last term of £4 in 
Eq. ([1]), for example, gives rise to a correction of the order 
arir2 to the bare gravitational constant. This affects the 
effective gravitational coupling Geff that appears in the 
equation of matter perturbations. Also it will be possible 
to constrain the Galileon models from the time variation 
of Gcff . Second, the study of spherically symmetric solu- 
tions in both weak and strong gravitational backgrounds 
can allow us to understand how the Vainshtein mecha- 
nism works in general. Third, it is possible to extend the 
Galileon Lagrangian ([Ij to the theory in which the field 
(f) is replaced by a general function /(^). We expect that 
such analyses will provide us deep insight on the possible 
modification of gravity and that it will shed new light on 
the nature of dark energy. 
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